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Abstract 



We discuss the relationship between the recurrence coefficients of orthogonal polynomials with respect to a semi- 
classical Laguerre weight and classical solutions of the fourth Painleve equation. We show that the coefficients in these 
recurrence relations can be expressed in terms of Wronskians of parabolic cylinder functions which arise in the description 
' of special function solutions of the fourth Painleve equation. 

1 Introduction 

In this paper we are concerned with the coefficients in the three-term recurrence relations for orthogonal polynomials with 
respect to the semi-classical Laguerre weight 

o 

CV u(x;t) = x x exp(-x 2 +tx), (1.1) 

with parameters A > — 1 and i e R, which has been recently studied by Boelen and van Assche [8] and Filipuk, van 
Assche and Zhang [23]. It is shown that these recurrence coefficients can be expressed in terms of Wronskians that arise 
. in the description of special function solutions of the fourth Painleve equation (Piv) 

" = i ( P) \ |« 3 + W + 2(, 2 - A)q + — , (1.2) 

dz 2q \ dz J 2 q 

where A and B are constants, which are expressed in terms of parabolic cylinder functions. 

The relationship between semi-classical orthogonal polynomials and integrable equations dates back to the work of 
Shohat [65] and later Freud [35]. However it was not until the work of Fokas, Its and Kapaev [27, 28] that these equations 
were identified as discrete Painleve equations. The relationship between semi-classical orthogonal polynomials and the 
(continuous) Painleve equations was demonstrated by Magnus [51, 52] who showed that the coefficients in the three-term 
recurrence relation for the Freud weight [35, 74] 

oj(x;t) = exp (-\x A ~ tx 2 ) , i£l, (1.3) 

with ieMa parameter, can be expressed in terms of solutions of Pry (1.2). 

A motivation for this work is the fact that recurrence coefficients of semi-classical orthogonal polynomials can often 
be expressed in terms of solutions of the Painleve equations. For example, recurrence coefficients are expressed in terms 
of solutions of Pn for semi-classical orthogonal polynomials with respect to the Airy weight 
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oj(x; t) — exp (ix 3 + tx) , x 3 < 



1 



with tela parameter [51]; in terms of solutions of P m for the perturbed Laguerre weight 

oj(x; t) — x a exp (—a; — t/x) , x € K + , 

with a > and t G R + parameters [12]; in terms of solutions of P v for the weights 

u(x; t) = (1 - x) a (l + xfe- tx , x e [-1, 1], 
w(a;;i) =a; a (l-a;)' 3 e-*/ x , a:e[0,l], 
w(a;;t) = x a \x + t) p eT x , x e M+, 

with a, /3 > and f e M+ parameters [3, 10, 34]; and in terms of solutions of Pyi for the generalized Jacobi weight 

u(x;t) =x a (l-xf(t-x)~>, xe[0,l], 

with a, P, 7 > and t e K+ parameters [20, 51]. 

Recurrence coefficients for orthogonal polynomials with respect to discontinuous weights which involve the Heaviside 
function %(x) have also been expressed in terms of solutions of Painleve equations [2, 14, 31, 33], while recurrence 
coefficients for orthogonal polynomials with respect to discrete weights have been expressed in terms of solutions of 
Painleve equations [7, 6, 18, 21, 22]. 

This paper is organized as follows: in §2, we review some properties of orthogonal polynomials; in §3, we review 
some properties of the fourth Painleve equation (1.2), including its Hamiltonian structure §3.1, Backlund and Schlesinger 
transformations §3.2 and special function solutions §3.3; in §4 we express the coefficients which arise in the three-term 
recurrence relation associated with orthogonal polynomials for the semi-classical Laguerre weight ( 1 . 1 ) in terms of Wron- 
skians that arise in the description of special function solutions of Pry (1.2); in §5 we derive asymptotic expansions for the 
recurrence coefficients; in §6 we discuss orthogonal polynomials with respect to the semi-classical Hermite weight 

oj(x;t) = \x\ x cxp(-x 2 + tx), x,teR, A > -1, 

which is an extension of the semi-classical Laguerre weight (1.1) to the whole real line, and show that the recurrence 
coefficients are also expressed in terms of Wronskians that arise in the description of special function solutions of Prv 
(1.2); and in §7 we discuss our results. 



2 Orthogonal polynomials 

Let P n (x), n € N, be the monic orthogonal polynomial of degree n in x with respect to a positive weight oj(x) on (a, b), a 
finite or infinite interval in R, such that 



J a 



P m (x)P n {x)uj(x)dx = h n 6 m ,n, h n >0, 



(2.1) 



where <5 TOi „ denotes the Kronekar delta. One of the most important properties of orthogonal polynomials is that they satisfy 
a three-term recurrence relationship of the form 

xP n {x) = P n+1 (x) + a n P n {x) + p n P n ^(x), (2.2) 

where the coefficients a n and f) n are given by the integrals 

1 f b 1 f b 

a n = — xP n (x)u(x)dx, f3 n = — / xP n ^ 1 (x)P r , 

"to Jo "n—l J a 

with P_i(x) = and Pq(x) — 1. These coefficients in the three-term recurrence relationship can also be expressed in 
terms of determinants whose coefficients are given in terms of the moments associated with the weight uj(x). Specifically, 
the coefficients a n and (3 n in the recurrence relation (2.2) are given by 



(x) uj(x) dx, 



(2.3) 



a r , 



+i 



+i 



An 

A„' 



A 



where A„ is the Hankel determinant 



A„ = det 



J j,k=0 



Mo Mi 

Ml M2 
Mn— 1 Mn 



ri+l^n-l 

Mn-1 
Mn 

M2n-2 



(2.4) 



n > 1, 



(2.5a) 
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with A = 1, A_i = 0, and A„ is the determinant 



A„ 



Mo Mi 
Mi M2 

Mn — 1 t^n 



fJ-n-2 

Mra-l Mn+1 
M2n-3 M2n-1 



n > 1, 



(2.5b) 



with A = and /Xfe, the fcth moment, is given by the integral 

rb 

Mfc 



= / X k Uj(x) 
J a 



dx. 



We remark that the Hankel determinant A n (2.5a) also has the integral representation 



^ rb rb n 

A„ = — / „• / || (xj - x k ) 2 dx 1 ... dx n , 



n > 1. 



(2.6) 



(2.7) 



The monic polynomial P„ (a;) can be uniquely expressed as the determinant 



Pn{x) = ±- 



Mo Mi 
Mi M2 



Mn— 1 Mn 
1 X 



Mn 
Mn+1 



M2n-1 



(2.8) 



and the normalisation constants as 



/in = 



-*n+l 
An 



h = Ai = ^o- 



For further information about orthogonal polynomials see, for example, [15, 43, 69]. 
Now suppose that the weight has the form 

w(x;t) — ljq(x) exp(xt), i£ [a, b], 



(2.9) 



(2.10) 



where t is a parameter, with finite moments for all (el, which is the case for the semi-classical Laguerre weight (1.1). If 
the weight has the form (2.10), which depends on the parameter t, then the orthogonal polynomials P n (x), the recurrence 
coefficients a n , (3 n given by (2.4), the determinants A„, A„ given by (2.5) and the moments ^ given by (2.6) are now 
functions of t. Specifically, in this case then 



J a 



fAk = I % uq(x) exp(xi) dx = 



u}q(x) exp(xt) dx ] = — -r~- 
/ dt 



Further the recurrence relation has the form 

xP n (x;t) = P n+1 (x;t) + a n (t)P n (x;t) + (3 n (t)P n - 1 (x;t), 
where we have explicitly indicated that the coefficients a n (t) and j3 n (t) depend on t. 



(2.11) 



Theorem 2.1. If the weight has the form (2.10), then the determinants A„(t) and A n (t) given by (2.5) can be written as 
the Wronskians 



A„(t)=W Mo 



dMo d ra Vq 

' dt ' ' ' ' ' J-"- 1 ' 



dt" 



A„(t) 



dA» 
dt ' 



(2.12) 



d fc Mo 



Proof. Since fi^ = — t~> tnen tne determinant A n (t) can be written in the form 

dt 



A„(i) 



Mo Mi • • • Mn-1 

Ml M2 ••• Mn 



Mn— 1 Mr. 



M2n-2 



= W mo 



d/^o d™ >o 



dt 



dt 



n-l / ' 
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as required, and the determinant A„(i), can be written in the form 

A„(t) 



Mo Mi 

Ml M2 



Mn-2 Mn 
Mra-1 Mn+1 



Mn— 1 Mn 



duo 



M2n-3 M2n-1 

d""Vo 



di 



d"-Vo d> 



»n — 1 



di ' 



dt ' ' d^ 
as required. 

The Hankel determinant A„ (i) satisfies the Toda equation, as shown in the following theorem. 
Theorem 2.2. The Hankel determinant A„ (t) given by (2.12) satisfies the Toda equation 

<P_ _ A n -i(t)A n+ i(t) 

dt 2 nW " A£(i) 

Proof. See, for example, Nakamira and Zhedanov [57, Proposition 1]; also [63, 67]. 



□ 



(2.13) 



□ 



Using Theorems 2. 1 and 2.2 we can express the recurrence coefficients a n (t) and f3 n (t) in terms of derivatives of the 
Hankel determinant A„ (t) and so obtain explicit expressions for these coefficients. 

Theorem 2.3. The coefficients a n (t) and f3 n (t) in the recurrence relation (2.11) associated with monic polynomials or- 
thogonal with respect to a weight of the form (2.10) are given by 



d A n+ i(t) 

a " ( * ) = d7 ln ^T 



(3 n (t) = -^lnA„(i), 



(2.14) 



with A n (t) is the Hankel determinant given by (2.12). 

Proof. By definition the coefficients a n (t) and (3 n (t) in the recurrence relation (2.1 1) are given by 

_ A n+1 (t) _ AJfl _ A n _ 1 (t)A n+1 (t) 

an(t) A n+1 (t) A n (ty Pn[t) ~ Al(t) 

where the determinants A„ and A„ are given by (2.5). Hence from (2.12) 

A n+l (t) A n (t) 1 dA„ +1 1 dA„ 



a n (t) 



and so 



as required. By definition 



and so from Theorem 2.2 we have 



A n+1 (t) A n (t) A„+i dt A„ dt 



d A n+ i(t) 

Q " w = d7 ln ^T 



(2.15) 



Pn(t) = 



A ra _!ft)A ra+1 ft) 



dt: 



■lnA„(t), 



as required. 

Additionally the coefficients a„ (i) and /?„ (t) in the recurrence relation (2. 1 1) satisfy a Toda system. 



(2.16) 
□ 



Theorem 2.4. The coefficients a n (t) and (3 n (t) in the recurrence relation (2.11) associated with a weight of the form 
(2.10) satisfy the Toda system 

= Pn+l - Pn, = Pn(a n ~ «n-l)- (2-17) 

Proof. See Ismail [43, §2.8, p. 41] and Moser [55]; see also [6] for further details and a direct proof in the case of the 
semi-classical Laguerre weight (1.1). □ 
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Suppose P n (x), for n € N, is a sequence of classical orthogonal polynomials (such as Hermite, Laguerre and Jacobi 
polynomials), then P n {x) is a solution of a second-order ordinary differential equation of the form 

d 2 P HP 

ff(2;) l^ i+r(x) ^ L = AnjP "' (2 ' 18) 

where a(x) is a monic polynomial with deg(er) < 2, t(x) is a polynomial with deg(r) = 1, and A„ is a real number 
which depends on the degree of the polynomial solution, see Bochner [5]. Equivalently, the weights of classical orthogonal 
polynomials satisfy a first-order ordinary differential equation, the Pearson equation 

^[a{x)w{x)] =t(x)lj(x), (2.19) 

with <j(x) and t(x) the same polynomials as in (2. 18), see, for example, [1, 5, 15]. However for semi-classical orthogonal 
polynomials, the weight function uo(x) satisfies the Pearson equation (2.19) with either deg(cr) > 2 or deg(r) > 1, see, for 
example, [41, 53]. For example, the Pearson equation (2.19) is satisfied for the weight (1.1) with 

a(x)=x, t(x) = -2x 2 + tx + A + 1, 

and so the weight (1.1) is indeed a semi-classical weight function. Filipuk, van Assche and Zhang [23] comment that 

"We note that for classical orthogonal polynomials (Hermite, Laguerre, Jacobi) one knows these recurrence 
coefficients explicitly in contrast to non-classical weights". 

In §4 we show that, in the case of the semi-classical Laguerre weight (1.1), the determinants A n (t) and A n (t) can be 
explicitly written as Wronskians which arise in the description of special function solutions of Piv (1-2) that are expressed 
in terms of parabolic cylinder functions D v (z) when A Z, or error functions erf (z) when A = n e Z. Consequently 
the recurrence coefficients a n (t) and f3 n (t) (2.4) associated with orthogonal polynomials for the semi-classical Laguerre 
weight (1.1) can also be explicitly written in terms of these Wronskians. 



3 Properties of the fourth Painleve equation 

The six Painleve equations (Pi-Pyi) were first discovered by Painleve, Gambier and their colleagues in an investigation of 
which second order ordinary differential equations of the form 

where F is rational in dq/dz and q and analytic in z, have the property that their solutions have no movable branch 
points. They showed that there were fifty canonical equations of the form (3.1) with this property, now known as the 
Painleve property. Further Painleve, Gambier and their colleagues showed that of these fifty equations, forty-four are 
either integrable in terms of previously known functions (such as elliptic functions or are equivalent to linear equations) or 
reducible to one of six new nonlinear ordinary differential equations, which define new transcendental functions, see Ince 
[42]. The Painleve equations can be thought of as nonlinear analogues of the classical special functions [17, 26, 39, 44, 71], 
and arise in a wide variety of applications, for example random matrices, cf. [30, 62]. 



3.1 Hamiltonian structure 

Each of the Painleve equations Pi-Pyi can be written as a Hamiltonian system 

dq = dUj dp = dHj 
dz dp ' dz dq ' 

for a suitable Hamiltonian function Hj(q,p,z) [45, 58, 60]. The function <r(z) = Hj(q,p,z) satisfies a second-order, 
second-degree ordinary differential equation, whose solution is expressible in terms of the solution of the associated Pain- 
leve equation [45, 59, 60]. 

The Hamiltonian associated with Prv (1.2) is 

H w (q,P, z; t?o, tfoo) = 2qp 2 - (<? 2 + 2zq + 2i9 )p + ^9, (3.3) 
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with ?? and ^oo parameters [45, 58, 59, 60], and so from (3.2) 

P- =4qp- q 2 -2zq - 20 o , (3.4a) 
dz 

^ = -2p 2 + 2qp + 2zp-$ 00 . (3.4b) 
dz 

Solving (3.4a) for p and substituting in (3.4b) yields 

dz 2 2q \dz J q 

which is Piv (1.2) with A = 1 - # + 2^ and 5 = -2t?Q. Analogously, solving (3.4b) for q and substituting in (3.4a) 
yields 

§ -5 (e) 2 + 6 " s - ^ + 2 <" - 2 "» + f • 

Then letting p = -±io yields Pi V (1.2) with A = -1 + 2i? - ^oo and B = -2i9^. 

An important property of the Hamiltonian, which is very useful in applications, is that it satisfies a second-order, 
second-degree ordinary differential equation. 

Theorem 3.1. Consider the function 

a(z; 0o, tfoo) - 2qp 2 - (q 2 + 2zq + 2i9 Q )p + ^g, (3.7) 
where q and p satisfy the system (3.4), then a satisfies the second-order, second-degree ordinary differential equation 

d 2 er\ 2 , ( da \ 2 da f da „ „ \ /da 



dz 2 ) -\ z te-°) +4 ^U +2 NU +2l?oo J =0 - (3 - 8) 

Conversely, if a is a solution of (3.8), then solutions of the Hamiltonian system (3.4) are given by 

_ a" - 2za' + 2a _ a" + 2za' - 2a , _ d 

q ~ 2(</ + 2# 00 ) ' P ~ 4 (</ + 2i? ) ' ~ d7" ( } 

Proo/ See Jimbo and Miwa [45] and Okamoto [58, 59, 60]. □ 
Remarks 3.2. 

1. Equation (3.8), which is often known as Srv (or the Prv cr-equation), is equivalent to equation SD-I.c in the classi- 
fication of second order, second-degree ordinary differential equations with the Painleve property by Cosgrove and 
Scoufis [19], an equation first derived and solved by Chazy [9]. 

2. Theorem 3.1 shows that solutions of equation (3.8) are in a one-to-one correspondence with solutions of the Hamil- 
tonian system (3.4), and so are in a one-to-one correspondence with solutions of Prv (1.2). 

3. Equation (3.8) also arises in various applications, for example random matrix theory [32, 33, 48, 70]. 

3.2 Backlund and Schlesinger transformations 

The Painleve equations Pn-Pyi possess Backlund transformations which relate one solution to another solution either of 
the same equation, with different values of the parameters, or another equation (see [17, 24, 39] and the references therein). 
An important application of the Backlund transformations is that they generate hierarchies of classical solutions of the 
Painleve equations, which are discussed in §3.3. 

Backlund transformations for Prv (1.2) are given as follows. 

Theorem 3.3. Let q = w(z; A , B ) and qf = w(z; Af,Bf), j = 1, 2, 3, 4 be solutions o/Prv (1-2) with 



Af = i(2-2A,±3v/-2£ ), Bf = -±(l + A ± ^-2B ) 2 , (3.10a) 

Af = -i(2 + 2A ±3v/-2B ), Bf = -±(l - A a ± ^-2B ) 2 , (3.10b) 

A± = §-±A)T W~ 2B o, B f =-\{l- Av±y-2B ) 2 , (3.10c) 

A i = -1 - h A o T W~ 2B o, Bf = -|(-1 -Ao± ^-2B Q ) 2 . (3.10d) 
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(3.11a) 
(3.11b) 
(3.11c) 



Then 

go-go - 2zg =F V-2B 

2go ^ 
g + gg + 2zg T V-2g 

2go ' 

2(l-A)T^V^27jQ-)go 
90 g ±V 3 2~B " + 2zgo + g 2 ' 
2(1 + A ±1^27^)g 

94 ~ 9 ° + „> =n / — OR" 9 7 n ^' ( } 

g T v-^-do - ^zgo - g 

valid when the denominators are non-zero, and where the upper signs or the lower signs are taken throughout each 
transformation. 

Proof. See Gromak [37, 38] and Lukashevich [50]; also [4, 39, 56]. □ 

A class of Backlund transformations for the Painleve equations is generated by so-called Schlesinger transformations 
of the associated isomonodromy problems. Fokas, Mugan and Ablowitz [29], deduced the following Schlesinger transfor- 
mations 7£i-7?.4 for Pry. 



7?: 


it 


7?: 


if 


It'- 


it 




it 



Tlx: qAz;A u B 1 ) = ^ '- — -A \ -, (3.12a) 

7? . a(z . A R ) (JL + - 4 - g 2 - g 2 (g + 2zf 
IZ 2 ■ q 2 (z; A 2 , B 2 ) = -, !— — r , (3.12b) 



(g' + V-2B)- + (4A + 4-2V- 


-25 


)g 2 


-q 2 (q + 2z) 2 


2g (g 2 + 2zg - q> - 


V- 


25) 




(g' - y/-2B)~ + [AA - 4 - 2V- 


-25 




-q 2 {q + 2zf 


2g (g 2 + 2zq + q' - 


V- 


25) 




(g'- y/-2B)~ - (4A + 4 + 2V- 


-25 


)1 2 


-q 2 (q + 2z) 2 


2g (g 2 + 2zq - g' + 


V- 


2B) 




(g' + V-2BJ" + (4^-4 + 2V- 


-25 


)l 2 


-q 2 (q + 2z) 2 



V 7 2g (g 2 + 2zg + g' + V^W) 

where g = q(z; A, B) and 

(A 1 ,B 1 )= (a + 1,-±(2-v^25) 2 ) , (A 2 ,5 2 ) = (A- 1,-1(2 + v / ^25) 2 ) , (3.12e) 

(A 3 ,5 3 ) = (a + 1,-\(2 + V^2B) 2 ) , (A 4 ,B 4 ) = (A- 1,-1(2- V^2B) 2 ) ■ (3.12f) 

Fokas, Mugan and Ablowitz [29] also defined the composite transformations 72-5 = IZ1R3 and 72-7 = TZ 2 TZ 4 given by 

(g'-g 2 -2zg) 2 +25 

72 5 : q 5 (z; A 5 , 5 5 = — ^ 1 io/^m ' < 3 " 13a ) 

2g {g - 1 — 2zq + 2 (A + 1)} 

(g' + g 2 +2zg) 2 + 25 

^ A ^ = - 2 q{ qlq" + 2zq-2(A-l } y ^ 

respectively, where 

(A 5 ,B 5 ) = (A + 2,B), (A 7 ,B 7 ) = (A-2,B). (3.13c) 
We remark that 1Z$ and IZ7 are the transformations 7+ and 71, respectively, given by Murata [56]. 

3.3 Special function solutions 

The Painleve equations Pu-Pyi possess hierarchies of solutions expressible in terms of classical special functions, for 
special values of the parameters through an associated Riccati equation, 

^ = / 2 (z)g 2 + / 1 (z)g + .f (z), (3.14) 

where f 2 (z), fi(z) and fo(z) are rational functions. Hierarchies of solutions, which are often referred to as "one -parameter 
solutions" (since they have one arbitrary constant), are generated from "seed solutions" derived from the Riccati equation 
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using the Backlund transformations given in §3.2. Furthermore, as for the rational solutions, these special function solutions 
are often expressed in the form of determinants. 

Solutions of Pii-Pyi are expressed in terms of special functions as follows (see [17, 39, 54], and the references therein): 
for Pn in terms of Airy functions Ai(z) and Bi(z); for Pni in terms of Bessel functions J v (z) and Y v (z)\ for Piv in 
terms of parabolic cylinder functions functions D v (z); for Py in terms of confluent hypergeometric functions \F\{a; c; z) 
(equivalently Kummer functions M(a, b, z) and U(a, b, z) or Whittaker functions M K ^(z) and W K ^{z)); and for Pyi in 
terms of hypergeometric functions 2-^1(0, b; c; z). Some classical orthogonal polynomials arise as particular cases of these 
special function solutions and thus yield rational solutions of the associated Painleve equations: for Pni and Py in terms of 
associated Laguerre polynomials (z); for Piy in terms of Hermite polynomials H n (z); and for Pyi in terms of Jacobi 
polynomials P^^\z). 

Special function solutions of Piy ( 1 .2) are expressed in in terms of parabolic cylinder functions. 
Theorem 3.4. Piy (1.2) has solutions expressible in terms of parabolic cylinder functions if and only if either 

B = -2(2n+ 1 + eAf, (3.15) 



or 

B = -2n 2 , (3.16) 

with neZ and e = ±1. 

Proof. See [36, 38, 39, 40, 49, 50]. □ 
For Piy (1.2) the associated Riccati equation is 

P- = e(q 2 + 2zq) + 2u, e 2 = l, (3.17) 

dz 

with Piv parameters A = —e(u + 1) and B = —2i> 2 . Letting w(z) = lmp„(z) in (3.17) yields 

dz 

^-^ - 2ez^p- + 2ev<p v = 0. (3.18) 
dz 2 dz 

The solution of this equation depends on whether v e Z or v £ Z, which we now summarize, 
(i) If v £ Z then equation (3.18) has solutions 

w (z . F ,\{CMV2z) + C 2 D u {-V2z)}e^{\z 2 ), if e = 1, 

VA ' ] \{C 1 D^ 1 (V2z) + C 2 D^ 1 (-V2z)}cxp(-±z 2 ), if e = -1, ^ * } 

with C\ and C 2 arbitrary constants, where D V (Q is the parabolic cylinder function which satisfies 

k¥Dv ={\t 2 -v-\)D v , (3.20) 



dc 2 V4S 

and the boundary condition 

A,(C)^rcxp(-iC 2 ), as C^+^- 
(ii) If v — then equation (3.18) has the solutions 

, , fG\+C 2 crt(z), if e = l, 
¥>o(z;e) = i (3.21) 
I 61 + 62 erfc(2), if e = — 1, 

with C\ and C 2 arbitrary constants, where erfc(z) is the complementary error function and erfi(z) is the imaginary 
error function, respectively defined by 

2 f°° 2 r z 

erfc(z) = -= / cxp(-i 2 ) dt, erfi(z) = -= / exp(t 2 ) dt. (3.22) 
V n Jz V n Jo 
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(iii) If v = to, for m > 1, then equation (3.18) has the solutions 



ip m (z;e) = < 



C\H m (z) + C 2 cxp(z 2 )-j- 7?r {crfi(z) cxp(-z 2 )} , 



C 1 (-i) m H m (iz) + C 2 cxp(-z 2 )^ {erfc(z)exp(z 2 )} 



if e = 1, 
if £ = -1, 



with Ci and C 2 arbitrary constants, where H m {z) is the Hermite polynomial defined by 

jm 

H m (z) = (-1)™ cxp(z 2 ) H?;r exp(-z 2 ). 



(3.23) 



(3.24) 



(iv) If ^ = —m, for to > 1, then equation (3.18) has the solutions 



(p- m (z; e) = < 



jm-l 



C 1 (-i) m - 1 ^ m _ 1 (iz)exp(z 2 ) + C 2 — 3T{erfc(z)exp( 2 2 )}, if £ = 1, 

, dz 

d m 

C l iJj" m _i(z)cxp(-z 2 ) + C 2 m-1 {erfi(z)exp(-z 2 )} , if £ = -1, 

dz 



(3.25) 



with C\ and C 2 arbitrary constants. 
If <p v (z; e) is a solution of (3.18), then the "seed solutions" of Piv (1.2) are given by 

q(z; —e(v + 1), -2v 2 ) = In <Pv(z; e), <?( 2; ; -e^ ~2(y + I) 2 ) = ~2z + e^- ]ntp„(z;e). 

Hierarchies of special function solutions can be generated from these solutions using the Backlund transformations given 
in §3.2. However there is an alternative approach. 

Determinantal representations of special function solutions for Prv (1.2) and Srv (3.8) are discussed in the following 
theorem. 

Theorem 3.5. Let T n ^(z; e) be given by 

T n , v {z;e) = W ^p v (z;e), ^-(z;e), . . . ,^-^f(z;s)j , n > 1, 



(3.26) 



with to, v (z; e) = 1, where <p v (z; e) is a solution of (3.18) and W(<pi,<p2, ■ ■ ■ , <Pn) is the Wronskian. Then for n > 0, 
special function solutions o/Prv (1-2) are given by 



e(-;4i,<i)=-2. + 4 ln ^ff: 

V / dz T ntI/ (z; 

Inju ( Z i ^n,vi Bn}Sj = £ ~T~ In ■ 



A z '-> £ ) 

(z;e) ' 



dz T n ^ +1 (z;e)' 

a [3] ( z . A [3] B [3] \ _ d m -r n ,„+i(z;e) 
v / dz T n +i,v(z;e) 

and special function solutions o/Srv (3.8) are given by 

Vnl( z ;$o,fioo) = -^lnT n ^(z;e)-2enz, 



AW V = s(2n - y), B$ v = -2(y + l) 2 , (3.27a) 
Ag} v = e(2v - n), B%} v = -2(n + l) 2 , (3.27b) 
4 3 l = s(n + v), 4% - -2(1/ - n + I) 2 , (3.27c) 



4 ] =e(v-n+l), tfW 



-en, 



)[2] _ 



= en, 



(3.28a) 
(3.28b) 



aW v (z-,# ,# oo ) = —lnT n , u (z-,s) + 2e(v-n + l)z, ^=-e(u + l), tfgj = -e(i/ - n+ 1), (3.28c) 



Proof. See Okamoto [60]; also Forrester and Witte [32]. 



□ 



4 Semi-classical Laguerre weight 

In this section we consider monic orthogonal polynomials P n (x;t), for n e N, with respect to the semi-classical Laguerre 
weight (1.1), where these polynomials satisfy the three-term recurrence relation (2.1 1), i.e. 



xP n (x;t) = P n+ i(x;t) + a n (t)P n (x;t) + fi n (t)P n -\(x;t), 
Boelen and van Assche [8, Theorem 1.1] prove the following theorem. 



(4.1) 
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Theorem 4.1. Let a n (t) and (3 n (t) be the coefficients in the recurrence relation (4.1) associated with the semi-classical 
Laguerre weight (1.1). Then the quantities 



V2 



t — 2a„ 

satisfy the discrete system 



2M = 2/3„-n-±A, (4.2) 



Vn + n + |A 1 / t 1 \ 

^n-l^n = 5 J-T^-, Vn + Vn+l = \ ~~ • (4-3) 

Boelen and van Assche [8] also show that the system (4.3) can be obtained from an asymmetric discrete Piv equation 
by a limiting process. However, from our point of view, it is more convenient to have the discrete system satisfied by a n 
and (3 n , which is given in the following Lemma. 

Lemma 4.2. The coefficients a n (t) and f3 n (t) in the recurrence relation (4.1 ) associated with the semi-classical Laguerre 
weight (1.1) satisfy the discrete system 

(2a n - t)(2a„_i -t) = - , (4.4a) 

Pn 

2/3 n + 2j3 n+l + a n (2a n - t) = 2n + A + 1. (4.4b) 
Proof. Substituting (4.2) into (4.3) yields the discrete system (4.4). □ 

Since the semi-classical Laguerre weight (1.1) has the form uio(x) exp(xi) and the moments are finite for all t 6 R, 
with t a parameter, then the coefficients a n (t) and f3 n (t) in the recurrence relation (4.1) satisfy the Toda system, recall 
Theorem 2.4. 

We are now in a position to prove the relationship between the coefficients a n (t) and /?„ (t) in the recurrence relation 
(4.1) associated with the semi-classical Laguerre weight (1.1) and solutions of Prv (1.2). 

Theorem 4.3. The coefficients a n (t) and j3 n (t) in the recurrence relation (4.1) associated with the semi-classical Laguerre 
weight (1.1) are given by 



a n (t) = \q n (z) + §t, (4.5a) 
A.W = -1% - k^n(z) + \n + \\ (4.5b) 



with z=\t, where q n {z) satisfies 



d 2 q n 1 / dq n \ , , 2 2 2A 



+ i9„+4z^ + 2(z 2 -2n-A-l) 9 „ , (4.6) 

dz 2q n \ dz J q n 

which is Piv (1-2), with parameters 

(A, B) = (2n + A + 1, -2A 2 ). (4.7) 

Proof. Solving the discrete system (4.4) for a n _i and (3 n +i yields 

i (2/3»-n)(2/?»-n-A) 
a "" 1 -^+ 2(2o n -i)/J„ 
Pn+i = ~Pn - |(2n + A + 1) - a n (a n - \t), 

and then substituting these into (2.17) gives 

= -a„(a„ - |t) - 2/?„ + i(2n + A + 1), (4.8a) 



Solving (4.8a) for /3„ yields 



d£„_, _i,, / o (2/3 n -n)(2/3 n -rc-A) 
dt 2 W " 2(2a n -i) 



(4.8b) 



(4.9) 
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and then substituting this into (4.8b) yields 

d 2 a n If da n 1 



(1 \ 2 

-^-h) + K - \ ta l + l( t2 - 4n - 2 - 2\)a n + \t(2n + A + 1) - 



dt 2 2a n -t\ dt V 2 4 4 4 4(2a„-t)' 

Making the transformation (4.5a) in this equation yields equation (4.6), which is Pry (1.2) with parameters given by (4.7). 
Finally making the transformation (4.5a) in (4.9) yields (4.5b), as required. □ 

Remarks 4.4. 

1. Filipuk, van Assche and Zhang [23], who considered orthonormal polynomials rather than monic orthogonal poly- 
nomials, proved the result (4.5a) for a n (t). However Filipuk, van Assche and Zhang [23] did not give an explicit 
expression for a n (t), which we do below. 

2. From Theorem 3.4 we see that the parameters (4.7) satisfy (3.15) with e = —1, and therefore satisfy the condition 
given in Theorem 3.4 for Prv to have solutions expressible in terms of parabolic cylinder functions. 

3. If q n is a solution of equation (4.6) then the solutions q n +i and q n -i are given by 
(q' n -ql-2zq n ) 2 - 4A 2 (q' n + q 2 % + 2zq n f - 4A 2 



<l " 1 1 2g„ (q' n - q 2 n - 2zq n + An + 2A + 4) ' 9n 1 2q n (q' n + q 2 n + 2zq n -An- 2A) ' 
where' = d/dz, which are special cases of the Schlesinger transformations 1Z§ (3.13a) and IZ7 (3.13b), respectively. 
4. From Theorem 3.5, we see that the parabolic cylinder function solutions of equation (4.6) are given by 

q n (z) = -2z+± ln^M (4.10) 



VW=W(^17 ^t], t ,a(z) = 1, (4.11) 



where 

T_ ,(z ) = W ( l/A. 

dz dz 
and t^\{z) satisfies 

^-2*^-2^ + 1)^=0, (4.12) 
dz 2 dz 

which is equation (3.18) with v = — A — 1 and e = 1. Equation (4.12) has general solution 

({C 1 D_ x _ 1 (V2z)+C 2 D- X _ 1 (-V2z)}cxp(±z 2 ), if A £ N, 

ipx(z) = I d m (4.13) 

[C 1 (-irH m (iz)exp(z 2 ) + C 2 — {erfc(z)exp(z 2 )}, if A = m e N, 

with Ci and C2 arbitrary constants, where D U {C) is the parabolic cylinder function, H m (() the Hermite polynomial 
(3.24), and erfc(z) the complementary error function (3.22). 

The system (4.8) satisfied by the recurrence coefficients a n (t) and P n (t) is equivalent to the Hamiltonian system (3.4) 
associated with Piv, as shown in the following Theorem. 

Theorem 4.5. The system (4.8) is equivalent to the Hamiltonian system (3.4) associated with Piv- 
Proof. If in the system (4.8) we make the transformation 

a n{t) = \q n (z) + \t, P n (t) = -\q n (z)p n (z) + \(n + \), z = \t, 
then q n (z) andp„(z) satisfy the system 

^ = 4q nPn -q\- 2zq n - 2A, (4. 14a) 

^ = -2p 2 n + 2 Pn q n + 2z Pn — n — X, (4. 14b) 

dz 

which is the system (3.4) with ^ = A and ?9oo = A + n. Conversely making the transformation 

. . . . , . 2f3 n {t) —n — X 

q n (z) = 2a n (t) - t, Pn{z) = — -- — — , t = 2z, 

2a n (t) - t 

in the system (4. 14) yields the system (4.8). □ 
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Our main objective is to obtain explicit expressions for the coefficients a n (t) and (3 n (t) in the recurrence relation (4. 1). 
First we derive an explicit expression for the moment ^o{t\ A). 

Theorem 4.6. For the semi-classical Laguerre weight (1.1), the moment fio(t; A) is given by 

r(A + l)exp(^) , y a ^ 



W>(*;A)=< J£ +1)/2 ' " ' (4-15) 

i\/^^{exp(it 2 ) [l + erf(|t)]}, if A = m e N, 

vwf/i D„ (£) f/ie parabolic cylinder function and erf (z) f/ze error function. Further fj, (t; A) satisfies the equation 

^-^-|(A + 1)mo = 0. (4.16) 
Proof. The parabolic cylinder function D„((), with !/^Z, has the integral representation [61, §12.5(i)] 

A,(C) = CX ^^ 2) a""" 1 exp(-i S 2 - C«) ds, R(iz) < 0. 
1 ( _!/ J Jo 

For the semi-classical Laguerre weight (1.1), the moment fio(t; A), with A N, is given by 

/■00 

/zo(i;A) = / x A exp(— x 2 + xt) Ax 
Jo 

= 2-( A + 1 )/ 2 y 0O S A cxp(-i S 2 + iV2t S ) ds 

r(A + l)exp(lt 2 ) / x 

= ^TT)75 V -2V2*J 

as required. IfmSN, then the parabolic cylinder function D_ m _i(Q is given by 

l»- m -i(0 = y|^rcxp(-IC 2 )^{ex P (K 2 )erfc(i^c)}, 

with crfc(z) the complementary error function [61, §12. 7(H)] . Since erfc(— z) = 1 + crf(z), then /J,o(t; m), with meN, 
is given by 

d m 

«,(t; m) = — {cxp (it 2 ) [1 + erf (it)] } , 

as required. Further, the parabolic cylinder function D V (Q) satisfies equation (3.20) and so from (4.15) it follows that the 
moment Ho(t; A) satisfies equation (4.16), as required. □ 

Corollary 4.7. Ifno(t; A) is given by (4.15) and <p„(z; e) by (3.19), then 

/io(t;A) = p_ A -i(§t;l), (4.17) 

with Ci = and C 2 = T(A + 1)/2( A+1 )/ 2 . 

Proof. The result is easily shown by comparing (4.15) and (3.19). □ 

Having obtained an explicit expression for ^ we can now derive explicit expressions for the Hankel determinant A„ (t) 
and the coefficients a n (t) and /3 n (t) in the recurrence relation (4.1). 

Theorem 4.8. The Hankel determinant A„ (t) is given by 



An(t)=W( W> ^,...,^ , (4.18) 



d^o d" Vo 
At At 

with no given by (4.15). 

Proof. This is an immediate consequence of Theorem 2.1. □ 
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Theorem 4.9. The coefficients a n (t) and f3 n (t) in the recurrence relation (4.1) associated with monic polynomials orthog- 
onal with respect to the semi-classical LagueRre weight (1.1) are given by 



d A n+ i(i) 
a n (t) = — In ■ 



n (t) = lnA„(t), 



(4.19) 



dt A n (t) ' 

where A n (i) is f/ze Hankel determinant given by (4.18), with fig given by (4.15). 

Proof. This is an immediate consequence of Theorems 2.1 and 2.3. □ 
Furthermore we can relate the Hankel determinant A„(t) given by (4.18) to the r-function r„ e) given by (3.26). 
Theorem 4.10. If A„(t) is given fey (4.75) anof r„^(z; e) fey (3.26), with 



r(A + l)cxp(lz 2 ) 
¥>-A-i(«) = n ao-iW9 



f/zen 



A„(t) 



2 (A+l)/2 

r„ 1) 



(-V2z) 



2«(n-l) 

Proof. The result is easily shown by comparing (4.18) and (3.26). 



z=t/2 



(4.20) 
□ 



Theorem 4.11. The function S n (t) = — In A n (t), with A n (t) given by (4.18), satisfies the second-order, second-degree 
equation 



d 2 S r , 
dt 2 



dS n 



^ dS n f^dSjt 
dt I dt 



dt 



Proof. Setting v = — A — 1 and e = 1 in (3.28c) gives 



cr(z;A,n + A) = — In T n -x-i (z; 1) - 2(n + A)z, 



(4.21) 



and so if (t) = — In A„ (t) then from (4.20) we see that 

a(z: : A, n + A) = 2S n (t) - (n + A)t, 
Making this transformation in Srv (3.8) with $ = \ and = n + A, i.e. 



z = it. 



(4.22) 



dV 

dz 2 



. dtr \ ^dcr / da 
\ dz / dz V dz 



+ 2A 



da 



+ 2n + 2A = 



yields (4.21), as required. 



□ 



Remark 4.12. Differentiating (4.21) and letting S n (t) = — In A„(t) yields the fourth-order, bi-linear equation 



A d 4 A» _ 4 _d 3 A„ dA„ 



dt 4 



dt 3 dt 



d 2 A r 
dt 2 



- (it 2 + An + 2A) 



A, 



d 2 A T , 
" dt 2 



dA T 
dt 



+ itA„^ + In(n + A)A 2 =0, 



as is easily verified. 

Theorem 4.13. Suppose ^ n ,\(z) is given by 



*„,*(*) = W Ua, 



dji\ 
dz 



in— 1 



dz T 



*o,a(*) = 1, 



where 



£-A-i(-\/2z)exp(iz 2 ), if A^N, 

^a(-z) = "( d m 

— |[l + crf(z)]cxp(z 2 )}, if A = me 



(4.23) 
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with -D;/(C) is parabolic cylinder function and erfc(z) the complementary error function (3.22). Then coefficients a n (t) 
and (3 n (t) in the recurrence relation (4.1) associated with the semi-classical Laguerre weight (1.1) are given by 

a n (t) = \q n (z) + \t, (4.24a) 

W!//z z =\t, where 

w/z/c/z satisfies P IV f^-2J, vvz'f/z parameters (A, B) = (2n + A + 1, — 2A 2 ). 

In Appendix 1 we give the first few recurrence coefficients for the semi-classical Laguerre weight (1.1) and the first 
few monic polynomials generated using the recurrence relation (4.1). 

5 Asymptotic expansions 

In this section we derive asymptotic expansions for the the moment [iq (i; A), see Lemma 5. 1 below, the Hankel determinant 
A„(i), see Lemma 5.2 below, and the recurrence coefficients a n (t) and /?„(<), see Lemma 5.3 below. 

Lemma 5.1. As f ^ oo, the moment Ho(t; A) has the asymptotic expansion 

*^*^wifti^iSiF (5 - i> 

Proof. Since the parabolic cylinder function -D„(C) has the asymptotic expansion 

Uv ^> r(-i,)C +1 PU4 J ^ n!(2C 2 )"' ^ ' 

with (/3)„ = T(/3 + n)/r(/3) the Pochhammer symbol, then 

r(A + l)exp(|t 2 ) / ^ x 

r(A + l)exp(lt 2 ) x/2^t A exp(^ 2 ) ~ (-A) 2n 

2 (A+l)/2 T(A + l)2 A /2 2^ n \ t 2n 

v ' n=0 

n=0 v ' 



as required, since 



. A , r(2n-A) W1 . „ „. T(A + 1) 

( - A)2 " = VaT = ( - 1) • ■ • (A - 2n + 1) = ^ J iy 



□ 



Lemma 5.2. As t — > oo, f/ze Hankel determinant A„ (i) /zas f/ze asymptotic expansion 

A n (t) = c„7r"/ 2 (Ii)" A exp ( ±rrf 2 ) j 1 + nA( *~ " } + O (fT 4 ) | , (5.2) 

with c n a constant, and S n (t) has the asymptotic expansion 

„ , . nA 2n\(n — A) _ , r N 

= y + — + ^ + O (t- 5 ) . (5.3) 
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Proof. To prove (5.2) we shall use Mathematical induction. Since A„ satisfies the Toda equation (2.13) then 

1 f A d 2 A» /dA„^- 
A n _! A " At 2 \ At 



A. + i = x^<lA^-(^) . (5.4) 



By definition A = 1 and from (5.1) 

A x = no = \pK {\tf exp (^ 2 ) 1 1 + + O (t- 4 ) | . (5.5) 

as t — > oo, and so (5.4) with n = 1 gives 



.^-(£)Vw»o^(K){. + ^ 



as £ — ^ oo. Assuming (5.2) then 



I di 2 \ dt / [ "2 n v 2 7 ^_p V2 .^ ; | j- ^ 2 
as i — > oo, and so 

— zb{^-(^)"} 

= ^L 7r (™+l)/2 ( l^(n+l)A exp | 1 (n + 1}i 2| 
^ c n— 1 

x { 1 + 2A(nA-n»-l) + Q ^ 4) j | x _ (n - 1)A(A - n + 1) + Q ^ j 

= c„ +l7 r("+ 1 )/ 2 (it)(«+ 1 )^ exp {!(„ + l)t 2 } {l + ( " ± 1)A(A = " = 1} + O (i~ 4 ) } 

as t — »■ oo, where c„ + i = inc 2 /c„_i, as required. Solving the recurrence relation 

c„+ic„_i = inc 2 , c = 1, ci = 1, 

gives 

^ n— 1 

C " = 2 n(n-l)/2 II kl 

Since 5„ = In A„ then the asymptotic expansion (5.3) is easily derived from (5.2). □ 
at ' " 

Lemma 5.3. As t — > oo, the recurrence coefficients a n (t) and (3 n {t) have the asymptotic expansions 

t A 2A(2n - A + 1) , 

M*) = g + 7 + — ^3 " + ° (* ) ' ( 5 - 6a > 

n , x n n ^ 6n\(n — A) _ , fiN 

=2^7^ 73 " + ° (* ) ■ ( 5 - 6b > 

Proof. By definition 

«nW = ^^% t 7^=S n+ i(t)-S„(t), n (t) = -^liiA n (t) <LS " 



and so 

as i — > 00, as required. Consequently 

lim a„(i) = ht, lim /3 n (i) = hn. 

t=oo z i=oo z 



□ 
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6 Semi-classical Hermite weight 

In this section we are concerned with the semi-classical Hermite weight 



oj(x;t) = |x| A exp(-x 2 +tx), i,tgK, A > — 1, (6.1) 

which is an extension of the semi-classical Laguerre weight (1.1) to the whole real line, where we have ensured that the 
weight is positive by using |x| A rather than x x . Monic orthogonal polynomials associated with the semi-classical Hermite 
weight (6.1) satisfy the recurrence relation 

xP n (x;t) = P„+i(x;t)+a„(f)P„(x;i)+/?„(t)P„_i(x;t), (6.2) 

and our interest is in obtaining explicit expressions for the coefficients a n (t) and /?„(t) in (6.2). 
First we evaluate the moment fj,o(t; A). 

Theorem 6.1. For the semi-classical Hermite weight (6.1), the moment fJ,o(t; A) is given by 

V^(-ii) 2m ^ 2m (iit)ex P (^ 2 ), if A = 2m, (6.3) 

_ j2m+l 

-^-2^+r{ erf (|*) exp(ii 2 )}, «f A = 2m + 1, 



po{t:X) = < 



df 



wif/i m € N, where D v (z) is the parabolic cylinder function, H n (z) is the Hermite polynomial and crf(z) is the error 
function. 

Proof. If A ^ N, then the moment ^o(t; A) is given by 

/oo />oo 
w(x;t)dx = / |x| A exp(— x 2 + tx) dx 
-OO J —CO 

/>OC P CO 

= / x x exp(— x 2 + tx) dx + / x x exp(— x 2 — tx) dx 
Jo Jo 

= r( * + 2 '^'V >-( - + <") 

as required. If A = 2m, with m 6 N, then 

/oo 
x 2m exp(-x 2 + tx) dx = ( - |i) m #2m(|tt) CX P (i* 2 ), (6-5) 
-OO 

as required, since the Hermite polynomial, H n (z), has the integral representation 

9 771 i'OC 

H m (z) = -i=\ (z + ix) m exp(-x 2 ) dx. 

V 71 " J -co 

Finally if A = 2m + 1, with m e N, then 

/oo 
x 2m \x\ cxp(-x 2 +tx) dx 
-co 

j2m / rco r-co 

= ^£2m [J x exp(—x 2 + tx) dx + J x exp(— x 2 — tx) dx 

/ roo rco \ 

= 2m+i ( / exp(— x 2 + tx) dx — / cxp(— x 2 — tx) dx) 

j2m+l 



— |0F {1 + erf(it)} exp (it 2 ) - {l - erf(|t)} exp (it 2 



as required, since 



_ j2m+l 

V^TFtMs') exp (it 2 )}, 



/•OO 

/ exp(-x 2 + tx)dx = iV7r {l + erf (it)} exp (it 2 
Jo 



□ 
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Next we obtain an explicit expression for the Hankel determinant A n (t). 
Theorem 6.2. The Hankel determinant A n (t) is given by 



where p,o(t; A) is given by (6.3). 

Proof. By definition the moment fj,k{t; A) is given by 

/oo 
a; fe |:r| A exp(— x 2 + tx) Ax 
-oo 

= ^(/°°N A ex P (-.x 2 +te)dx) dV ° 



and so we obtain 



A n (t) = det fj, j+ k(t) 



W [no 



dt k 

du d n ~Vo 



j,k=o ' ' V" 1 " dt ' ■ ' ■ ' dt 71 - 1 
as required. □ 

Finally we obtain explicit expressions for the coefficients a n (t) and (3 n (t). 

Theorem 6.3. The coefficients a n (t) and (3 n (t) in the recurrence relation (6.2) associated with monic polynomials orthog- 
onal with respect to the semi-classical Hermite weight (6.1) are given by 

a «W = ^ ln % ± 77T' p n (t) = ^lnA n (t), (6.7) 
at A„(tj dt 

where A n (t) is the Hankel determinant given by (6.6), with fio{t; A) given by (6.3). 

Proof. This is an immediate consequence of Theorem 2.3. □ 
Theorem 6.4. Suppose ^> n ,\{z) is given by 



where 



'{D_A-i(v^«) +D_a-i(- V2z)}exp(iz 2 ), if A £ N, 
V>a(z) = l #2mO)exp(> 2 ), if A = 2m, m G N, 

— 2^TT { erf (' z ) ex P( 2;2 )} ' if A = 2m+1, rneN, 

■ dz 

with D V (Q is parabolic cylinder function, H m (C) the Hermite polynomial (3.24), and erfc(z) the complementary error 
function (3.22). Then coefficients a n (t) and f3 n (t) in the recurrence relation (6.2) associated with the semi-classical 
Hermite weight (6.1) are given by 

a n (t) = \q n (z) + ft, (6.8a) 
Pn(t) = - \ql(z) - \zq n {z) + \\+ in, (6.8b) 

with z = \t, where 

w/zzc/z satisfies Pry (7.2), vvz'f/z parameters given by (A, B) = (2n + A + 1, — 2A 2 ). 

In Appendix 2 we give the first few recurrence coefficients for the semi-classical Hermite weight (6.1), in the case when 
A = 2 (so the recurrence coefficients are rational functions of t), and the first few monic polynomials generated using the 
recurrence relation (6.2). 
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7 Discussion 



In this paper we have studied semi-classical Laguerre polynomials which are orthogonal polynomials that satisfy three-term 
recurrence relations whose coefficients depend on a parameter. We have shown that the coefficients in these recurrence 
relations can be expressed in terms of Wronskians of parabolic cylinder functions. These Wronskians also arise in the 
description of special function solutions of the fourth Painleve equation and the second-order, second-degree equation 
satisfied by the associated Hamiltonian function. Further we have shown similar results hold for semi-classical Hermite 
polynomials. The link between the semi-classical orthogonal polynomials and the special function solutions of the Painleve 
equations is the moment for the associated weight which enables the Hankel determinant to be written as a Wronskian. In 
our opinion, this illustrates the increasing significance of the Painleve equations in the field of orthogonal polynomials and 
special functions. 
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Appendix 1. Recurrence coefficients and polynomials for the semi-classical Laguerre weight 

For the semi-classical Laguerre weight the first few recurrence coefficients are given by 



£-A-i(-±V2t) 



a {t) = \t- \ 2 - J iX = tf„(i), 



ai(t) = |i-*„(*) *" (/) 



2*2( t ) _ _ A - L 



, . , 2A + 4 VJt) 



~>2 



_ 2[(A + l)t 2 + 4(A + 2)(2A + 3)]ttg(t) - (A + l)t[t 2 + 2(4A + 9)]tt„(t) — (A + l) 2 [t 2 + 8(A + 2)] 
2t[2t*3(i) - (t 2 - 4A - 6)*2(t) - 3 (A + l)**„(t) - 2(A + l) 2 ] 

{t) = -^l(t) + \W u (t) + \(\ + l) 1 

= 2t*l(t) - (t 2 - 4A - 6)* 2 (t) - 3(A + l)t*„(i) - 2(A + l) 2 



2[n{t)- 1 1 w u {t)-\(\+i)Y 

and the first few monic orthogonal polynomials are given by 

Pi{x;t) =x-^ v 

p . ™l - (t 2 + 2)tt„ - (A + l)t 2(A + 2)* 2 - (A + - (A + 1) 



P 3 {x;t) = x 3 - 



2[* 2 - \t* v - i(A + 1)] 2[* 2 - |t*„ - |(A + 1)] 

4(< 2 + 2A + A)^l - 2t{t 2 - A - l)* 2 - (A + l)(5t 2 + 4A + 6)tf „ - 3(A + ift ' 
2 [2t*jj - (t 2 - 4A - 6)*2 - 3(A + l)t* v - 2(A + l) 2 ] 



J 2i(t 2 + 2A + 4)*3 _ [(4 + 4 (2A + 5)(A + 2)] * 2 - 2(A + l)t(t 2 - A - 5)tf„ - (A + l) 2 (t 2 - 4A - 12) 1 
+ \ 4[2i*3 - (t 2 - 4A - 6)* 2 - 3(A + l)ttf„ - 2(A + l) 2 ] ) X 

2 [(A + l)t 2 + 4(A + 2) 2 ] - (A + l)t{t 2 + 2A + 8)^ 2 - 2(A + l) 2 (t 2 + 2A + 5)tt„ - (A + ift 
+ 4[2t*3 - (t 2 - 4A - 6)* 2 _ 3 (A + - 2(A + l) 2 ] 
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Appendix 2. Recurrence coefficients and polynomials for the semi-classical Hermite weight 

For the semi-classical Hermite weight x 2 exp(— x 2 + tx) the first few recurrence coefficients are given by 

a {t) = 

it 3 2t 



ai( * ) = ' t+ ?+12 

n _ i 6i(t 4 + 12 - 4i 2 ) 4t 3 
a 2 (ij - 2* + ;6 _ 6i 4 + 36t 2 + 72 - 44 + 12' 

n _i 8i 3 (t 4 + 60- 12t 2 ) 6t(t 4 + 12 - it 2 ) 

. , , 10i(t 8 + 216i 4 + 720 - 24i 6 - 480i 2 ) 8t 3 (t 4 + 60 - 12< 2 ) 



2 t w - 30t 8 + 360t 6 - 1200i 4 + 3600i 2 + 7200 t 8 - Wt 6 + 120t 4 + 720 

3 1^+8 _ An + 6 — i — _ , OAnn\ 1 f\+(+8 J- 91 fi/ 4 



a 6 (t) 

&(*) = 1 - 

&(*) = §- 

/3 4 (i) = 2- 



! 12t 3 (t 8 - 40i 6 + 600t 4 - 3360i 2 + 8400) 10i(t 8 + 216i 4 + 720 - 24i 6 - 480t 2 ) 

2 * + t 12 - A8t 10 + 900t 8 - 6720i 6 + 25200t 4 + 100800 ~ t 10 - 30t 8 + 360t 6 - 1200i 4 + 3600t 2 + 7200' 



1 2ft 2 -2) 

2 (t 2 + 2) 2 ' 
_ 4f 2 (t 2 -6)ft 2 + 6) 

(t 4 + 12) 2 

3 6(i 4 -12t 2 + 12)(i 6 + 6t 4 + 36< 2 -72) 



(i 6 - 6t 4 + 36t 2 + 72) 2 
8t 2 (t 4 - 20t 2 + 60) (t 8 + 72t 4 - 2160) 
ft 8 - 16i 6 + 120t 4 + 720) 2 ' 
10(i 6 - 30t 4 + 180i 2 - 120) (t 12 - 12i 10 + 180i 8 - 480i 6 - 3600i 4 - 43200i 2 + 43200) 
(t 10 - 30t 8 + 360t 6 - 1200i 4 + 3600t 2 + 7200) 2 ' 



and the first few monic orthogonal polynomials are given by 

t(t 2 +6) 



Pi(x;t) = 
P 2 {x;t) = x 2 - 



2(t 2 + 2)' 

t(t 4 + At 2 + 12) t 6 + 6t 4 + 36i 2 - 72 

t^Tn x + 4(t 4 + 12) ' 



, » , M(t 6 - 2t 4 + 20t 2 + 120) 3(t 8 + 40i 4 - 240) t(t 8 + 72t 4 - 2160) 



2(i 6 - 6t 4 + 36t 2 + 72) 4(< 6 - 6i 4 + 36t 2 + 72) 8{t 6 - 6t 4 + 36t 2 + 72) : 
4 2t(t 8 - 12t 6 + 72t 4 + 240i 2 + 720) 3 3{t w - 10t 8 + 80t 6 + 1200< 2 - 2400) 2 
4 (x;t)-x t8 _ — 6 + i2Qt4 - x + 2 (i 8 - 16t 6 + 120i 4 + 720) X 

t(t 10 - 10t 8 + 120i 6 - 240i 4 - 1200i 2 - 7200) 



2(t 8 - 16t 6 + 120t 4 + 720) X 
t i2 _ m io + 18Qt 8 _ 480 ^6 _ 36 oot 4 - 43200t 2 + 43200 
+ 16(t 8 - 16i 6 + 120t 4 + 720) ' 



P 5 (x;t) =x 5 - 



.4 



5t(10080 + t w + 264t 6 + 1680i 2 - 26t 8 - 336i 4 ) 
2ft 10 - 30t 8 + 360i 6 - 1200i 4 + 3600i 2 + 7200) 
5(t 12 - 24t 10 + 252t 8 - 672t 6 + 5040t 4 - 20160) 
2(t 10 - 30i 8 + 360i 6 - 1200t 4 + 3600i 2 + 7200) 
5t(t 12 - 24t w + 300t 8 - 1440t 6 + 5040i 4 - 100800) 
4(t 10 - 30i 8 + 360t 6 - 1200t 4 + 3600i 2 + 7200) 



•x 2 



5{t u - 26t 12 + 396t 10 - 2520i 8 + 5040t 6 - 50400i 4 - 100800i 2 + 201600) 
+ 16(t 10 - 30i 8 + 360t 6 - 1200t 4 + 3600t 2 + 7200) X 

t(t 14 - 30i 12 + 540t 10 - 4200t 8 + 10800t 6 - 151200i 4 - 504000i 2 + 3024000) 
32(t 10 - 30i 8 + 360i 6 - 1200t 4 + 3600i 2 + 7200) ' 
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